Scattering Properties of Paramagnetic Ground States in the 
Three-Dimensional Random-Field Ising Model 

Gaurav P. Shrivastav^, Siddharth Krishnamoorthy^, Varsha Banerjee^ and Sanjay Puri^ 
^School of Physical Sciences, Jawaharlal Nehru University, New Delhi - 110067, India. 
^Department of Physics, Indian Institute of Technology, 
Hauz Khas, New Delhi - 110016, India. 

Abstract 

We study the ground-state (T = 0) morphologies in the d = 3 random-field Ising model (RFIM) 
using a computationally efficient graph-cut method. We focus on paramagnetic states which arise 
for disorder strengths A > Ac, where Ac is the critical disorder strength at T = 0. These 
paramagnetic states consist of correlated "domains" of up and down spins which are separated 
by rough, fractal interfaces. They show novel scattering properties with a cusp singularity in the 
correlation function at short distances. 



PACS numbers: 64.60.De - Statistical mechanics of model systems: Ising model, Monte Carlo techniques, 
etc.; 68.35.Rh - Phase transitions and critical phenomena; 75.60.Ch - Domain walls and domain structure 



Spin systems with quenched disorder have challenged physicists for several decades. The 
competing nature of the interactions creates difficulties in studying them analytically and 
computationally. As a result, the properties of phases and phase transitions in disordered 
systems remain controversial. The random-field Ising model (RFIM) is an archetypal exam- 
ple of a system with quenched disorder and is described by the Hamiltonian l|, I2I: 

N 

E = ffiffj - ^ /liCTj, (Ji = ±1. (1) 

(ij) i=l 

Here, J > is the strength of the exchange interaction between nearest-neighbor spins. The 
variables {hi} are random fields, usually drawn from a Gaussian distribution whose standard 
deviation A is a measure of disorder. The phase diagram of the RFIM has been the subject 
of much discussion. In the 2-dimensional case {d = 2), there is no long-range order in the 
presence of disorder, no matter how small. However, in (i = 3, there is a small region of 
(T, A)-values where the equilibrium phase is ferromagnetic 0, 4|. Let us focus on the case 
with zero temperature (T = 0). In that case, the system exhibits a phase transition from a 
ferromagnetic phase (for A < Ac) to the paramagnetic phase (for A > Ac). The nature of 



this transition has received considerab 
RFIM is due to Middleton and Fisher 



An important study of the d = 3 



e attention 

^ . They studied a wide range of physical properties 
and convincingly demonstrated that there is a second-order phase transition at A = Ac. 

At T = 0, all the information about the system is encoded in the ground-state. Further, 
according to the zero-temperature fixed point hypothesis, transitions at T = and T 7^ are 
in the same universality class jo], [lo|. Therefore, a study of the ground-state morphology 
is important in understanding the RFIM phase diagram in c? = 3. A typical method of 
accessing the ground-state is via Monte Carlo (MC) evolution (e.g.. Metropolis [11], Simu- 
lated annealing etc.) from an arbitrary initial condition. However, MC approaches for 
disordered systems suffer from several drawbacks. First, the competition between exchange 
interactions and the random field introduces deep valleys in the free-energy landscape. These 
metastable states trap the evolving system and impede the relaxation to the ground-state. 
The system then opts for a local minimum, which can be far removed from the global min- 
imum, and may not refiect any of its properties. Further, as the MC techniques involve ~ 
order (1) spin-fiip at a time, the possibility of escape from a local minimum to the global 
minimum is small. Second, MC methods suffer from a non-polynomial (NP) divergence of 
computation time with system size. Thus, it is computationally very demanding to reach 



the global minimum for large systems with disorder. Consequently we still do not have a 
complete understanding of the nature of the ground-state. 

To address this problem, several optimization techniques based on "max-flow/min-cut" 



or "graph-cuts" have 
of binary variables 



Deen developed for a wide class of energy functions (or Hamiltonians) 
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15| . The basic approach in a graph-cut method (GCM) is to construct 



a specialized graph for the energy function to be minimized such that the minimum cut on 
the graph also minimizes the energy. The cut enables simultaneous relabeling of several 
spin variables or nodes. As a consequence, an exponentially large portion of the phase space 
can be sampled in a single move, thereby facilitating a quick search for a global minimum 
or a "good-quality" local minimum. Typically, the search time in these procedures has 
a polynomial dependence on the system size. An important class of energy functions are 
those which are (a) quadratic and (b) satisfy a "regularity" condition. In that case, the 
max-flow/min-cut technique actually yields the global minimum or exact ground state of 
the energy function in polynomial time IgI-IsI. The Hamiltonian of the RFIM specified in 
Eq. ([T]) belongs to this class [19|. We are therefore assured of reaching the exact ground 
state of the RFIM if energy minimization is via graph-cuts. 

The literature on combinatorial optimization provides many graph-cut algorithms with 
different polynomial complexity times. Some of the standard approaches include the Ford- 
Fulkerson (FF) method of augmenting paths I3!], the Goldberg- Tarj an (GT) push-relabel 
method [jjI, and the more recent Boykov-Kolmogorov (BK) method jl5|. A benchmarking 
of the above algorithms on a number of typical graphs has revealed that the BK method 
works several times faster than any other. While the FF and GT algorithms exhibit an 
A^^- dependence on the system size A^, the BK method is linear in A^ [15]. 

In this paper, we use the BK method to study the ground-state (T=0) morphologies 
in the d = 3 RFIM. This enables us to access exact ground states for substantially larger 
system sizes than in previous studies js]. We need these large sizes to obtain smooth data 
for statistical properties of the morphology, e.g., correlation function, structure factor, etc. 
We focus on the scattering properties of the domain structure in the paramagnetic state, i.e., 
for field strengths A > A^. This domain morphology has several non-trivial features, which 
we highlight in this paper. The correlation function C (r. A) (= {(XiCrj) — {ai){aj) with r = 
1?^ — 7^1) is a scaling function of r/^ (A), where the correlation length C, diverges as A — t- A^. 
At small values of r/^, C (r. A) exhibits a cusp singularity characterized by the roughness 
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exponent a: C{r,A) ~ 1 — A (r/^)" + ■ ■ -. This singularity has important consequences 
for the high-momentum behavior ("tail") of the structure factor. A similar cusp has been 



reported earlier in the context of fluctuation- dominated phase separation [20|, |2l|, and is a 
consequence of soft and ragged interfaces separating equilibrium phases. In the paramagnetic 
phase of the RFIM, there are no coexisting equilibrium phases. Nevertheless, there exist 
correlated domains of size ~ C, which are enriched in up or down spins. The scattering 
properties of these domain boundaries are analogous to those of fractal interfaces. We also 
provide accurate estimates of the critical point Ac and the correlation length exponent u 
calculated from the ground-state morphologies. 

Before presenting our results, we discuss the graph-cut approach, which has many po- 
tential applications for energy minimization in complex spin systems. This method can be 
applied to energy functions of the form: 

Ei{s.}) = E V,,is,,s,) + Y,D.i^^)- (2) 

The label Si of site i & S can take a value or 1, and the sites are related to one another by 
a well-defined neighborhood JV. The function Di measures the cost of assigning the label Sj 
to the site i, and Vij{si,Sj) measures the penalty (or cost) of assigning labels Si and Sj to 
adjacent sites i and j. 

The starting point in a GCM is to construct a specialized graph for the energy function 
E such that the minimum cut on the graph yields minimization of the energy. A graph Q is 
an ordered pair of disjoint sets (V, £), where V is the set of vertices and £ is the set of edges. 
An edge ij joining vertices i and j is assigned a weight Vij. A cut C is a partition of the 
vertices V into two sets TZ and Q. Any edge ij G S with i E TZ and j G Q (or vice-versa) is 
a cut edge. The cost of the cut is defined to be the sum of the weights of the edges crossing 
the cut. The minimum-cut problem is to find the cut with the smallest cost. 

The energy function E must satisfy the regularity condition for it to be graph- 
representable. Regularity is defined by the inequality Vij{0, 0) + V^j(l, 1) < Vij{l, 0) + Vij{0, 1) 



171]. The spin variables (cij = ±1) in Eq. ([T]) can be transformed into occupation-number 
variables (rij = 0, 1) through the transformation = (1 + cri)/2. Then, neglecting constant 
terms, 



N 



E{{ni}) = -AJ^niUj - 2'^{hi - qJ)ni, = 0, 1, (3) 

{ij) i=l 



where q denotes the number of nearest neighbors of a lattice site. It is straightforward to 
check that the interaction term of Eq. ([3]) satisfies the regularity condition. Thus, the energy 
function E of Eq. ([T]) is graph-representable and can be minimized using a GCM to yield 
the exact ground state. Each iteration of the GCM finds an optimal subset of nodes with a 
fixed label s, (= or 1) that gives the largest decrease in energy. This computation is done 
via graph-cuts on the specialized graph representing the energy function E. The algorithm 
repeatedly cycles through the labels Si until the global minimum is reached. 

Our simulations of the T = RFIM have been performed on d = 3 lattices of size 
{L < 256), with periodic boundary conditions applied in all directions. All the statistical 
data presented here, unless otherwise specified, is for a cubic system with L = 256. The 
initial configuration of the lattice is chosen to be a random mix of CTj = ±1, corresponding 
to the paramagnetic state at A = oo with ^ = 0. The results have been averaged over 100 
sets of {hi} for each value of A. Our studies indicate that the GCM has a 99% overlap with 
the ground-state in the first iteration itself, provided the disorder strength is not too close 
to the critical value Ac. (We do observe "critical slowing down" in the GCM as A — )■ Ac, 
but the phenomenon is much milder than in conventional MC methods.) We also find that 
the average energy per spin in the ground-state is an order of magnitude less than that 
obtained by the Metropolis algorithm. As mentioned earlier, the MC evolution invariably 
gets trapped in high-energy metastable states. 

The ground-state morphology in the paramagnetic state has the following features. As 
the disorder strength is reduced from A = oo, there is emergence of correlated regions or 
domains of size ^, enriched in either up or down spins (see snapshots in Fig. [T]for A = 2.4, 
2.6, 2.8). These regions grow in size with ^ — )■ oo as A — > A^. [A similar divergence of 
^ is seen for A = 0, T — )■ T+. However, in that case, a typical MC snapshot (see Fig. [2]) 
shows that the domain morphology is not as compact or well-defined as in Fig. [Tl F ig. [2] 
corresponds to T = 4.515, which is very close to the critical point: = 4.5103 22|.] By 
computing the Binder cumulant 23| as a function of A for different system sizes, we estimate 
Ac (T = 0) ~ 2.278 ± 0.002. For A < Ac, the ground-state morphology consists of a single 
large domain of ordered spins (up or down) with small clusters of oppositely-directed spins. 
The fraction of oppositely-directed spins decreases as A — )■ 0. 

What are the quantitative properties of the domain morphologies in Fig. [Tp The standard 
probe for quantifying these patterns is the correlation function C (r. A). The correlation 
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length ^(A) is defined as the distance over which C (r, A) decays to (say) 0.2 x maximum 
value. In Fig. [3t^a) we plot ^{5, L) vs. 5 [where 5 = ( A — Ac)/ Ac] for system sizes L ranging 
from 16 to 256. The correlation length diverges as ^ ~ 5^'^ when 5 but this is limited 

by the lattice size L. The presence of these finite-size effects can be used to estimate the 
critical exponent v. The finite-size scaling ansatz ^(5, L) = 5~^f{L5'^) results in the data 
collapse seen in Fig. |3]^b), yielding v ~ 1.308 ± 0.005. This is consistent with the earlier 
results of z/ = 1.4 ± 0.2 (Rieger and Young j2^) and v = 1.37 ± 0.09 (Middleton and Fisher 

1)- 

If the system is characterized by a single length scale, the morphology of the domains 
does not change with A, apart from a scale factor. In that case, the correlation function 
exhibits scaling: C (r. A) = g{r/C,) 25j. This is verified in Fig. HI where we plot C (r. A) 
vs. r/C, for different disorder amplitudes A > Ac. The data collapse for different values of 
A is excellent, confirming that the morphologies are scale-invariant. 

Next, we turn our attention to the central theme of this paper, viz., the scattering prop- 
erties of the domain morphology in Fig. [H The scattering of a plane wave by a rough surface 
can yield useful information about the texture of the surface Thus, small-angle scat- 

tering experiments (using X-rays, neutrons, etc.) can be used to probe the nature of domain 
walls separating the components of an inhomogeneous system. These experiments yield 
the structure factor S{k,A), which is the Fourier transform of the correlation function. 
Experimentalists are interested in the large-fc (tail) behavior of S{k,A), which is deter- 
mined by the small-r (short-distance) behavior of C (r. A). In the inset of Fig. HI we plot 
1 — C (r. A) vs. r/C, on a log-log scale. The small-r behavior shows a distinct cusp singularity: 

C (r. A) ~ 1 — A {r / C,)°' -\ with a ^ 0.5. This holds over more than a decade in r/^-values. 

A similar cusp has been reported earlier also in the context of fluctuation- dominated phase- 



separation (FDPS) 20|, |2l|. The cusp exponent a is identical to the roughness exponent 



of the domain boundaries. [Notice that the value we obtain for the paramagnetic phase 
(«para — 0.5) differs considerably from the roughness exponent in the ferromagnetic phase. 
The Middleton-Fisher value for the latter exponent is aferro = 0.66 ±0.03, which is consistent 



with the theoretical result afcrm = 2/3 29|, |30|.] The corresponding interfaces are self-affine 



fractals with df = d — a. Therefore, in our present study, dj ~ 2.5, which is consistent with 
percolation clusters in the strong-disorder regime of the d = 3 RFIM by Seppala 



studies o 
et al. 



3l| and Ji and Robbins [32|. 
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In Fig. [5l we show a schematic of a domain of size ^, with an interface of width w. There 
is also a microscopic length scale a = 1, due to the underlying discreteness of the lattice. The 
corresponding C (r, A) would show corrections to scaling, characterized by the parameter 
w/C,- Systems characterized by a cusp singularity exhibit very rough interfaces with w ~ ^. 
A novel feature of the present work is that we observe this scattering phenomenology in 
the paramagnetic phase of the RFIM, i.e., in the absence of interfaces between coexisting 
equilibrium phases. The lower frames in Fig. [T]show cross-sections of the d = 3 snapshots in 
the upper frames. Notice that the domain structure is fuzzy and subject to large fluctuations, 
and the boundaries are ill-defined. 

At larger values of a; = r/^, the correlation function is well-approximated as 

C (r. A) ^ 1 - A (r/0" - B (r/0 + ■ ■ (4) 
The linear decay in Eq. is characteristic of scattering from sharp interfaces in inhomoge- 



neous systems, and is termed the Porod law 33| . With reference to the schematic in Fig. 0, 



the correlation function C(r, A) exhibits (a) no systematic structure for r ~ a; (b) interfacial 
structure or cusp singularity ioi w ^ r ^ a; (c) Porod decay for ,^ ^ r ^ t;;. 

The short- distance cusp singularity in C(r, A) has important implications for the struc- 
;ure factor S {k, A). The scattered intensity now decays with an asymptotic power-law form 



26|428| 



S {k, A) ~ i (e^)"^''^"^ + B i^k)-^''^'^ , (5) 

valid for k ^ a^^. For k ~ a~^, the structure factor becomes flat, corresponding to the 
absence of structure at microscopic scales. The dominant large-fc behavior in Eq. ([5]) is 
S{k) ~ with cross-over momentum fcc ~ ^ • 

In Fig. El we plot the structure factor for the RFIM with A > Ac on a log-log scale. 
Our data is consistent with the scaling form in Eq. ([5]). There is a cross-over from a Porod 
regime [with S {k, A) ~ k^^'^'^^^~\ at intermediate values of k to an asymptotic cusp regime 
[with 5* {k, A) ~ k^^'^^°'\ a ~ 0.5] . The inset of Fig. [6] shows the behavior of the cross-over 
momentum, which scales as kc ~ We conjecture that the crossover is a generic feature 
in the RFIM as a consequence of interfacial roughening caused by quenched disorder 3J] . In 
this context, we consider Refs. [35|,l36| where the authors studied domain growth in the c? = 3 
RFIM. They focused on the nonequilibrium evolution of the system after a quench from the 



paramagnetic phase (A = oo) to the ferromagnetic phase (A < Ac). Refs. [35|, l36| observe 



that the scahng functions of the RFIM and the pure Ising system are identical, thereby 



exhibiting super-universality 



37H39| or irrelevance of quenched randomness. However, a 



careful observation of the scaled correlation data (Fig. 2 in 35|) for small r/^ reveals clear 
deviations from the pure system and so from the Porod law. [Of course, we expect to 
recover the Porod law in the limit w/C, — 0. This is possible in the domain growth problem 
as ^(t) — > oo as t — )■ oo.] 

We conclude this paper with a summary and discussion of our results. We have used 
a computationally efficient graph-cut method (GCM) to study the ground-state (T = 0) 
properties of the RFIM in the paramagnetic state. The Boykov-Kolmogorov GCM used by 
us provides access to the ground-state morphology of large systems. We characterize this 
morphology using correlation functions and structure factors, which contain information av- 
eraged over all domains and interfaces. The correlation function C (r. A) is characterized by 
a universal scaling function for different disorder amplitudes. There are no perceptible cor- 
rections to scaling for different values of A, suggesting that the interface thickness w scales 
with the correlation length ^. At short distances, C (r. A) shows a cusp singularity reminis- 



cent of that seen in fluctuation- dominated phase separation 20|, l2l|]. This is associated with 



scattering off rough, fractal interfaces. The corresponding structure factor S {k, A) shows a 
crossover from a Porod regime at intermediate k values, to an asymptotic cusp regime. These 
properties should be universal for disordered systems, which are often characterized by rough 
interfaces. We believe that our results will motivate further analytical and numerical studies 
of this problem. 
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A= 2.4 



A = 2.6 



A = 2.8 





FIG. 1: Ground-state morphologies of the RFIM obtained using the a-expansion GCM, for disorder 
strengths A = 2.4, 2.6 and 2.8. The snapshots in the top frames correspond to a 64^ lattice with 
periodic boundary conditions in all directions. Regions with up spins and down spins are marked 
black and grey. The domains shrink in size and interfaces roughen with increasing disorder, as is 
evident from the cross-sections (taken at z = 32) in the bottom frames. 
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FIG. 2: Equilibrium morphology of the disorder-free Ising paramagnet for T = 4.515 
(Tc ~ 4.5103 for A = 0). This snapshot is obtained from Glauber MC simulations of a 64^ lattice 
with periodic boundary conditions. The system is evolved from an arbitrary initial condition to its 

equilibrium state, where the morphology is invariant with time. The above snapshot corresponds 
to a cross-section (at z = 32) of a 64^ lattice. 
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FIG. 3: (a) Plot of correlation length ^{S,L) vs. S, where S = {A — Af.) / A^. We present data for 
cubic lattices of size L = 16, 32, 64, 128, 256 - denoted by the specified symbols. The data sets were 
averaged over 100 random-field configurations. The correlation length is defined as the distance 
over which the correlation function falls to 0.2 x maximum value, (b) Data collapse resulting 
from the finite-size scaling ansatz i{6^L) = d~'^f{L6'^), yielding the correlation length exponent 
u ~ 1.308. 
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FIG. 4: Scaled correlation function [C (r, A) vs. r/^] for specified disorder strengths. The numer- 
ical data has been averaged over 100 random- field configurations for a lattice of size 256^. The 
inset shows the small-r/^ behavior on a log- log scale to highlight the cusp singularity. The slope 
of the solid line yields the cusp exponent a ~ 0.5. 
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FIG. 5: Schematic of a domain of size ^. The characteristic interface thickness is w, and the 
microscopic lattice spacing a = 1. 
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FIG. 6: Structure factor [S'(fc, A) vs. A:], corresponding to the correlation functions in Fig. 4. The 
solid lines denotes a Porod regime [S* {k, A) --^ /c"^] at intermediate A;-vahics, which crosses over 
to an asymptotic cusp regime [S {k, A) ~ fc"^-^] . The inset shows the behavior of the cross-over 
momentum kc vs. ^"^ for several values of disorder. 
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